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We propose and study a new class of of superconducting detectors which are sensitive to Cl(meV) 
electron recoils from dark matter-electron scattering. Such devices could detect dark matter as light 
as the warm dark matter limit, mx > 1 keV. We compute the rate of dark matter scattering off 
of free electrons in a (superconducting) metal, including the relevant Pauli blocking factors. We 
demonstrate that classes of dark matter consistent with terrestrial and cosmological/astrophysical 
constraints could be detected by such detectors with a moderate size exposure. 


Introduction. The search for the identity of dark 
matter (DM) is in an exciting and rapidly developing 
era. Theories of weakly interacting massive particles 
(WIMPs) for DM, being predictive and testable, have 
been the primary focus of both theory and experiment for 
the last thirty years. Strong constraints from direct de¬ 
tection experiments, such as XenonlOO [T], LUX [5] and 
SuperCDMS [3], along with the absence of new physics 
signals from the LHC, have, however, been painting such 
models as increasingly constrained and tuned. Further, 
because the energy threshold of direct detection experi¬ 
ments searching for WIMPs is typically 1-10 keV, these 
experiments lose sensitivity to DM particles with mass 
below 10 GeV. At the same time, DM candidates with 
low masses are theoretically well-motivated: asymmetric 
dark matter HIS] and strongly interacting massive par¬ 
ticles [5] are examples in which the natural mass scale of 
the DM sits beneath the ^ 10 GeV scale. 

A new frontier for massive DM thus opens for 1 keV < 
mx < 10 GeV, with the lower bound set approximately 
by warm dark matter constraints, e.g. from phase space 
packing [3 [S] or the Lyman-a forest [3]. For elas¬ 
tic scattering processes, the deposited energy is Ed — 

/(2me,N), where q ^ PrVx is the momentum trans¬ 
fer with vx ^ 10“^ the incoming DM velocity, /i^ the 
reduced mass of the system and is the mass of 

the target electron or nucleus N. Thus for 100 MeV 
DM, an eV of energy is deposited for scattering off a nu¬ 
cleus. Inelastic processes, such as electron ionization or 
excitation above a band gap, may occur when the DM 
kinetic energy exceeds the binding energy. Utilizing a 
semi-conducting crystal such as germanium, with a band 
gap of 0.7 eV, implies potential sensitivity to DM as light 
as C>(MeV) [101 [H]- SuperCDMS is already working to 
lower its threshold to 300 eV [3] , constraining 1 GeV mass 
DM. 

To go well below this, as low as the warm DM limit at 
C>(keV), requires a different kind of technology; in this 
case one must be able to access electron recoil energies 
as low as O(meV). The purpose of this letter is to inves¬ 
tigate a proof of principle experiment to search for DM 


down to the warm DM limit. Devices utilizing supercon¬ 
ductors, we will show, are ideal for this purpose, as they 
can be sensitive to extremely small energy depositions. 
In fact, in cold metals, the limit on the sensitivity of the 
experiment to low energy DM recoils is set by the ability 
to control the noise rather than by an inherent energy 
gap in the detector. 

The targets we discuss are metals, with the DM in¬ 
teracting with free electrons in the Fermi sea. The DM 
scattering rate is limited by Pauli blocking for electrons 
locked deep in the sea, yielding a suppression factor of 
order the energy transfer over the Fermi energy; the sup¬ 
pression is, e.g., of order ~ 10“^ for a DM-electron scat¬ 
tering with meV energy deposition in a typical metal such 
as aluminum. As we will show, DM models satisfying all 
astrophysical and terrestrial constraints are detectable 
despite the Pauli blocking effect, extending the concep¬ 
tual reach of the detection method down to DM masses 
of C>(keV). 

Detection with Superconductors. The chal¬ 
lenge in designing a detector to observe DM with low 
energy deposits is to achieve a large target mass, while 
keeping noise low. Detection of small energy deposi¬ 
tions is by now well-established; superconductors, with 
a meV superconducting-gap, have sensitivity to energies 
at this scale. Transition edge sensors (TES) and Mi¬ 
crowave Kinetic Inductance Devices (MKIDs) have been 
utilized to detect microwaves and x-rays with sub-meV 
to keV energies in astrophysical applications. For exam¬ 
ple, TESs with sensitivity to energy depositions not very 
far from our range of interest already exist: Refs. HMD 
have demonstrated noise equivalent power in the range 
~ 10“^® — 10 “^°W/-\/Hz. This translates to a sensitiv¬ 
ity of ^ 50 — 300 meV of energy over a read-out time of 
~ 10 ms. Thus current technology could already start 
probing new regions of parameter space, though not yet 
at the C>(meV) level of sensitivity required for probing 
down to keV dark matter. Since the energy resolution 
scales with VT^V, with T the heat bath temperature 
and V the TES volume, the required improvement could 
be made by lowering T (from, for instance, 100 mK to 
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10 mK) and further decreasing the heat capacitance of 
the TES by reducing the volume. 

The TES and MKID, however, have very low masses— 
an MKID is typically a nano-gram in weight, while TESs 
are approximately 50 microns on a side and a fraction 
of a micron thick. As a result, they do not make good 
detectors themselves. Their masses cannot simply be in¬ 
creased since this would decrease their sensitivity to small 
deposits of energy. An alternative is then to use the TES 
or the MKID merely as heat sensors which register small 
deposits of energy from a much larger target mass, an 
‘absorber.’ 

Eor the absorber, we choose a supeconductor; a super¬ 
conductor features an energy gap which controls the ther¬ 
mal noise in the absorber. As a DM particle hits a free 
electron in the Fermi sea of the absorber, the recoiling 
electron will deposit an 0 ( 1 ) fraction of its energy into 
breaking Cooper pairs, creating quasiparticles in the su¬ 
perconductor. These quasiparticles random walk in the 
superconductor until the energy stored in them can be 
collected. Two possibilities for the collection are that 
the quasiparticles (I) re-combine and create an athermal 
phonon or (II) are absorbed on collection fins on the sur¬ 
face of the absorber. 

In the former case, the athermal phonon may break 
Cooper pairs in the MKID, leading to an observable 
change in the kinetic inductance. In the latter case, the 
quasiparticles may reach a collection fin on the surface of 
the absorber. The fins should have a lower gap than the 
absorber, both to control noise and to facilitate collection 
of energy into the fins. The collection fins are connected 
to the TES which registers the heat. The quasiparticle 
lifetimes are sufficiently long and their velocities suffi¬ 
ciently high that even if the collection fin area on the 
absorber is small, the quasiparticles ricochet sufficiently 
many times that they are very efficiently channeled from 
the absorber into the collection fins and on to the TES. 
Aluminum is an example of an ultra-pure metal that 
makes for a good absorber: with quasiparticle lifetimes of 
order a millisecond [15] and velocities of order the Fermi 
velocity vp ^ 10“^c, its gap of ^ 0.3 meV pairs well 
with gapless gold collection fins. We note that the scat¬ 
tering length in the absorber sets the upper bound on 
its unit size—of order ^ 5 mm in ultra-pure aluminum— 
such that many small absorbers must be multiplexed for 
large exposure. 

In either case, the MKID or the TES is acting as a 
calorimeter for the energy deposited in the absorber. The 
underlying design principle sketched here is of concen¬ 
tration: one seeks to store the deposited energy non- 
thermally, whether through quasiparticles or athermal 
phonons, and then concentrate them through a collection 
mechanism onto the MKID or TES. This process must 
happen fairly rapidly, on the timescale of a millisecond. 

Our purpose here is not to advocate for a particular 
experimental design, but rather simply to outline how. 


through improvements to existing technology, sensitivity 
to extremely light DM utilizing superconductors may be 
feasible. (Other techniques, such as the use of superfluid 
helium m, hold promise as well.) The remainder of this 
letter focuses on the reach of such an experiment into the 
parameter space of light dark matter. 

Rates and Backgrounds. Detection via TESs (or 
MKIDs) operates by DM scattering off of free electrons 
in a metal. In a superconductor, the free electrons are 
bound into Cooper pairs, which typically have ~ meV 
(or less) binding energy. Once the energy in the scat¬ 
tering exceeds this superconducting gap, however, the 
scattering rate is computed by the interaction with free 
electrons, times a coherence factor. This factor is 0(1) 
for energies just above threshold, and goes to unity for 
energies above the gap, see e.g. Ref. m- In the setups 
we consider, the gap is below the noise-limited energy res¬ 
olution, and the coherence factor can be neglected. The 
electrons are then described by a Fermi-Dirac distribu¬ 
tion at low temperature. The typical Fermi energy Ep of 
these electrons is p\/(2me) 10 eV, with pp ~ 50 keV 

in a typical metal such as aluminum. Scattering with 
a target electron buried in the Fermi sea can break the 
Cooper pair if the energy transferred in the scattering 
is enough to pull an electron out of the sea and above 
the gap. As a result, with kinetic energy of the incoming 
DM approximately mxv\ ^ meV— keV for keV to GeV 
DM, Pauli blocking is important for the DM scattering 
rate. We follow the discussion in [18] to compute the 
rate correctly, factoring in the Pauli blocking effect. We 
denote the 4-momentum of DM initial and final states 
by Pi and P 3 , the initial and final states of the electron 
by P 2 and P 4 , and the momentum transfer q = (Ed, q). 
The scattering rate can be estimated via 
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where Ep is the deposited energy, (|A4p) is the squared 
scattering matrix element summed and averaged over 
spin, and fi(Ei) = [1 -|- exp()]~^ is the Fermi- 
Dirac distribution of the electrons at temperature T. 
S(Ed, |q|) characterizes the Pauli blocking effects, and 
in the limit of T —)■ 0 , ^(EDjIql) reduces to a simple 
Heaviside theta function, with amplitude mlEp/(Tr\q\). 
We perform the integral numerically in order to capture 
the entire kinematic range properly. The total rate (per 
unit mass per unit time) is then 


Ed 
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dEp 


J dvxfMB(vx)ED 


d(neO'Urel) 1 PX 
dEp p mx 


( 2 ) 


Here p is the mass density of the detector material, and 
px = 0.3 GeV/cm3 the DM mass density. We take the 
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velocity distribution of the DM /mb(^’x) to be a modified 
Maxwell Boltzmann with rms velocity vq = 220 km/sec, 
and cut-off at the escape velocity Vesc = 500 km/sec. 
Since the typical Fermi velocity of a metal is vp = 
0(10^) km/sec » Vesc, frei — vf- The Pauli blocking ef¬ 
fect provides a suppression factor of order Ejj/Ep, which 
we confirm numerically. An irreducible background is ex¬ 
pected to come from electron-neutrino scattering, which, 
due to the low energy deposition in the detector, will be 
dominated by pp neutrinos uniiin]. We find that the solar 
neutrino background is many orders of magnitude below 
the signals we consider, and is hence omitted from fur¬ 
ther discussion. We have also checked that backgrounds 
from Compton scatters (at levels already achieved in an 
experiment such as CDMS) are not significant. 

In what follows we assume that the DM X interacts 
with electrons via exchange of a mediator (j). The gen¬ 
eralization of light DM models will be addressed in fu¬ 
ture work 1^; we seek only to demonstrate proof of 
principle here. The scattering cross section between, 
e.g., Dirac DM and free electrons is given by Uscatter = 
167raeQ;xMex/(^| + where ai = gt is the 

coupling of (/) to i with i = e,X, peX the reduced mass of 
the electron-DM system, and q the momentum transfer 
in the process. This cross-section is related to the ma¬ 
trix element in Eq. Q via Uscatter = 

We define two related reference cross sections udd, cor¬ 
responding to the light and heavy mediator regimes: 


- light 
^DD 

- heavy 
^DD 


167raeQ!x 2 

t^eX ■ 

^ref 

167raeQ!x 2 

4 ^^eX ■ 

ml 


9ref = PeXVx , 


(3) 


where vx = 10“^. The transition between these regimes 
is set by how large the mediator mass is in comparison 
to the momentum transfer. The reference momentum 
transfer g^ef above is chosen for convenience as a typical 
momentum exchange. Note however that for a light me¬ 
diator, the direct detection cross section is determined by 
the minimal momentum transfer in the process, which is 
controlled by the energy threshold of the detector. 

To establish a notion of the expected number of events, 
in Fig. we present the differential rate per kg-year as 
a function of deposited energy for several benchmark 
points described in the next section. When the medi¬ 
ator is effectively massless, namely lighter than the mo¬ 
mentum transfer in the scattering, the rate is peaked at 
energies near the detector threshold due to the 1/g^ en¬ 
hancement of the cross section. In contrast, for mas¬ 
sive mediators, the rate is peaked at higher recoil en¬ 
ergies. The reason for the latter behavior is that as the 
recoil energy increases, more electrons can be pulled from 
deeper in the Fermi sea, resulting in an increased rate. 
The mass of the mediator determines the scattering dis¬ 
tribution in phase space, but does not control the size 



FIG. 1: Signal rates per kg-year, for several benchmark points 
of {m^, mx, ax,ge) = (10 /reV, 10 keV, 5 x 10“^^^, 3 x 10“®) 
[solid green], (10 geV, 100 MeV, 5x10“®, 3x 10“^®) [dashed 
green], (1 MeV, 10 keV,0.1,3 x 10“®) [solid red], and 
(100 MeV, 100 MeV, 0.1,3 x 10“®) [dashed blue[. We use 
the Fermi energy of aluminum, Ep = 11.7 eV. 


of the available phase space. A cutoff in the differen¬ 
tial rate is evident for both light and heavy mediators, 
and depends on the DM mass. For heavier DM (dashed 
curves), the maximum energy deposition is determined 
by £1)5“ = \me{{vF + 2nesc)^ — v'p)- When the DM 
is lighter (solid curves), the cutoff is determined by the 
kinetic energy of the DM, namely by /rex'Cesc/2- 

Results. In Fig. we show the 95% expected sensi¬ 
tivity reach after one kg-year exposure, corresponding to 
the cross section required to obtain 3.6 signal events 
The left (right) panel corresponds to the light (heavy) 
mediator regime, where we plot ^ func¬ 

tion of mx ■ The black solid [dashed] curve in both panels 
corresponds to a sensitivity to measured energies between 
1 meV—1 eV [10 meV—10 eV]. For light mediators, the 
scattering rate is sensitive to the lowest energy deposi¬ 
tions, resulting in a large improvement in reach when 
the detector threshold is decreased. For massive media¬ 
tors, the differential rate peaks towards larger energies, 
though with a lower threshold there is more sensitivity 
to lighter particles. 

For a sense of the size of the cross-sections in Eq. (|^, 
we divide our discussion into light mediator and heavy 
mediator regimes. We begin with a light mediator cj), 
which for illustration purposes we take to be a scalar. 
In the left panel of EigJ^ we plot the direct detection 
cross section dQ§* [Eq. (1^ for several benchmark points 
labeled I-III, shown in solid colored curves. As is evi¬ 
dent, large direct detection cross sections can be obtained 
even for extremely small couplings due to the large en¬ 
hancement factor in Eq. (§, that scales like four powers 
of the inverse of the momentum transfer in the detection 
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Massive mediator 



FIG. 2: Left: Direct detection cross section for light dark matter scattering off electrons, for several benchmarks of light 
mediators. These are I: ax = 10“^®,Oe = II: ax — a^ = and III: ax ~ 10“^®,Oe = These depicted 

parameters obey bounds from self-interactions and decoupling at recombination for < eV, though stellar emission may 
place strong constraints for scalar mediators; see text for details. Right: Direct detection cross section between light dark 
matter and electrons, for several benchmarks of heavy mediators. These are A: = 1 MeV, Qe = 10“®e, ax ~ 0.1; B: 

= 10 MeV, ge ~ 10“®e, ax = 0.1; and C: — 100 MeV, ge ~ 10“^e, ax = 0.1. These depicted parameters obey 

all terrestrial and astrophysical constraints, though sub-MeV DM interacting with SM through a massive mediator may be 
strongly constrained by BBN; see text for details. The XenonlO electron-ionization data bounds m are plotted in thin dashed 
gray. In both panels, the black solid (dashed) curve depicts the sensitivity reach of the proposed superconducting detectors, 
for a detector sensitivity to recoil energies between 1 meV—1 eV (10 meV—10 eV), with a kg-year of exposure. For comparison, 
the gray dot-dashed curve depicts the expected sensitivity utilizing electron ionization in a germanium target as obtained in 
Ref. [lO]. 


process when the mediator is light. The presented bench¬ 
mark points all obey DM self-interaction bounds 
and also ensure that the DM remains out of kinetic equi¬ 
librium with the baryons up through the time of recom¬ 
bination [2H] for < eV. Stellar constraints are model- 
dependent (for example, whether a scalar or kinetically 
mixed vector mediator), and hence have not been fac¬ 
tored in here; we note that for a kinetically mixed hid¬ 
den photon, the strength of stellar constraints is lifted 
for the couplings shown in the plot since the combina¬ 
tion ~ gx9e or m 0 ( 7 e is then bounded [SD] rather 
than just g^. Also note that the reach curves do not 
include any medium-dependent mediator mass, as this 
is model-dependent. For example, in a metal, a kinet¬ 
ically mixed vector mediator would experience a large 
in-medium mass; such a mass becomes small in an in¬ 
sulating superfluid absorber like helium. We detail the 
medium- and model-dependence in a longer paper Ell- 

Moving to heavy mediators, we focus on > MeV. 
A plethora of constraints exists in the literature for this 
mass range, see e.g. EIH311 in the context of kinetically 
mixed hidden photons. In the right panel of Fig. we 
select several benchmark points, labeled A-C, that sur¬ 
vive all terrestrial {e.g. beam dump) and stellar cooling 
constraints, and plot the resulting direct detection cross 
section of Eq. (|^, Large couplings to electrons 

ge > 10 “® are possible despite stellar constraints due 


to trapping effects, and beam dump constraints may be 
evaded by decaying to additional particles in the dark 
sector. These statements hold regardless of the vec¬ 
tor/scalar nature of the heavy mediator. However, for 
values of ax and ge as large as these benchmark points, 
DM and/or the mediator will be brought into thermal 
equilibrium with the SM plasma. The chief constraint on 
these models is thus BBN and Planck limits on the num¬ 
ber of relativistic species in equilibrium (see e.g. |35)). 
The Planck constraints can be evaded; for instance cou¬ 
pling to y/e through the time that the DM becomes 
non-relativistic will act to reduce the effective number 
of neutrinos at CMB epoch. On the other hand, dur¬ 
ing BBN, the helium fraction constrains the Hubble pa¬ 
rameter, which is sensitive to all thermalized degrees of 
freedom. DM must then be either a real scalar or heav¬ 
ier than a few hundred keV in such simple models [35) . 
It follows that part of the depicted curves of benchmarks 
A-C in the low-mass region may not be viable; a detailed 
study of the viable parameter space is underway [2T] . For 
completeness, we show the XenonlO electron-ionization 
bounds [30] in the thin gray dashed curve. (The XenonlO 
bounds on light mediators are not depicted in the left 
panel of Fig. as they are orders of magnitude weaker 
than the parameter space shown.) 

For comparison, we show the expected sensitivity using 
electron-ionization techniques with a germanium target 
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as obtained in Ref. HU], translating their result into ctdd 
of Eq. These results are depicted by the dot-dashed 
gray curves in Fig. for both the light (left panel) and 
heavy (right panel) mediator cases. For heavy media¬ 
tors and mx larger than a few hundred keV, our de¬ 
tection method is less sensitive than the projected one 
using germanium, while for lighter mx, where electron 
ionization methods lose sensitivity, the superconducting 
devices win. (Indeed, this comparison between the de¬ 
tection methods is our main aim in presenting the right 
panel of Fig. [^) In contrast, light mediators highlight the 
strength of our proposed detectors. For DM masses above 
several hundred keV, superconducting detectors can out¬ 
perform electron ionization techniques by several orders 
of magnitude. For dark matter below the MeV scale, the 
proposed superconducting detectors are uniquely staged 
to detect super light sub-MeV viable models of dark mat¬ 
ter. 

In summary, we have proposed a new class of detec¬ 
tors that utilize superconductors to detect electron recoils 
from thermal DM as light as a keV. Given some improve¬ 
ment over current technology, such detectors may have 
sufficiently low noise rates to be sensitive to the required 
energy scale of meV electron recoils. We have computed 
the DM scattering rates, taking into account Pauli block¬ 
ing, and have shown that viable models may be detected. 
We hope this proof of concept encourages the experi¬ 
mental community to pursue research and development 
towards the feasibility of such devices, probing detection 
of DM down to the keV scale. We leave for future work 
the extended study of broader classes of DM models that 
may be detectible with these devices. 
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